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Abstract: We study the power- and bi-spectrum of vacuum fluctuations in a hyperbolic 
section of de Sitter space, comparing two states of physical interest: the Bunch-Davies and 
hyperbolic vacuum. We introduce a one-parameter family of de Sitter hyperbolic sections 
and their natural vacua, and identify a limit in which it reduces to the planar section and 
the corresponding Bunch-Davies vacuum state. Selecting the Bunch-Davies vacuum for a 
massless scalar field implies a mixed reduced density matrix in a hyperbolic section of de 
Sitter space. We stress that in the Bunch-Davies state the hyperbolic de Sitter re-point 
correlation functions have to match the planar de Sitter re-point correlation functions. The 
expressions for the planar and hyperbolic Bunch-Davies correlation functions only appear 
different because of the transformation from planar to hyperbolic coordinates. Initial state 
induced deviations from the standard inflationary predictions are instead obtained by con¬ 
sidering the pure hyperbolic vacuum, as we verify explicitly by computing the power- and 
bi-spectrum. For the bi-spectrum in the hyperbolic vacuum we find that the corrections 
as compared to the standard Bunch-Davies result are not enhanced in specific momentum 
configurations and strongly suppressed for momenta large compared to the hyperbolic cur¬ 
vature scale. We close with some final remarks, in particular regarding the implications of 
these results for more realistic inflationary bubble scenarios. 
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1 Introduction 

Cosmological observations point to a primordial universe that can be effectively described 
by an approximate de Sitter phase. The details of this inflationary phase are to a large 
extent still unknown, but the most recent Planck data does rule out a large fraction of 
parameter space, giving us some hints about the underlying physics [1]. From a theoret¬ 
ical point of view it is probably fair to say that inflation is poorly understood. Indeed, 
few convincing theoretical constraints on the inflationary parameter space exist that would 
identify natural, UV consistent, models. Trying to embed inflation into a fundamental 
description like string theory is notoriously difficult, but the last decade has seen consid¬ 
erable phenomenological progress in that direction (for a thorough discussion see [2] and 
references therein). The absence of a guiding principle which is able to rule out a significant 
fraction of inflationary models, combined with the attractive features of eternal inflation. 
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has fueled the idea that perhaps our inflationary universe is just one realization in a huge 
landscape of bubble universes that are continuously being produced as a consequence of a 
stochastically varying scalar field during a phase of eternal inflation [3, 4]. 

This exotic possibility makes the general prediction that the spatial sections in our 
universe should be hyperbolic on the largest scales that start probing the boundary of 
the bubble. Clearly this general prediction is hard, if not impossible, to verify because the 
primordial inflationary expansion typically redshifts the negative curvature scale far beyond 
the observable universe [5], although in the inflationary landscape relatively short phases 
of slow-roll inflation might be preferred [6], which could lead to observable consequences 
in the CMB temperature correlations at low multipoles [7]. More particular predictions 
include a bubble universe that might, under fortuitous conditions, provide an explanation 
for the low power anomaly at low I due to the steepening of the slope right after penetration 
of the barrier [8]. Another potential consequence of observational interest is that bubble 
universes can in principle collide [9-11], which could leave definite non-isotropic signatures 
in the Cosmic Microwave Background sky. Unfortunately the chance for a collision to have 
taken place in our past is small and model-dependent, and therefore not seeing this effect 
will not be able to rule out a multiverse origin [12]. 

Recently, another rather generic consequence of a multiverse origin has been explored 
[13-15]. In the context of an inflationary landscape one would expect the initial vacuum 
state for quantum fluctuations in a single inflationary bubble to be entangled with the 
rest of the eternally inflating universe, leading to a mixed state inside the bubble. Since 
the Cosmic Microwave Background temperature anisotropies (as well as the large scale 
structure distribution) is probing the statistics of these inflationary quantum fluctuations, 
one could imagine uncovering evidence in favor of a mixed initial state that would support 
the idea that our universe originated from false vacuum decay. This idea warrants a careful 
study of the actual (observational) potential to constrain departures away from a standard 
pure initial state for inflationary quantum fluctuations and how these departures relate to 
the global vacuum state of eternal inflation. 

In this work we present a first step, triggered by some recent work in this direction 
[13], in clarifying the connection between the vacuum state of the false, eternally inflat¬ 
ing, vacuum and potential departures from the standard Bunch-Davies vacuum state in 
a hyperbolic bubble. As a starting point we identify a limit that connects hyperbolic 
coordinates and its naturally associated vacuum to the de Sitter invariant Bunch-Davies 
state on planar sections. Following up on older work in the context of open inflation, after 
selecting the global Bunch-Davies vacuum we then use the mixed reduced density matrix 
defined in a single hyperbolic coordinate patch of de Sitter space to explicitly show that 
the statistics of inflationary quantum fluctuations are indistinguishable from the standard 
planar Bunch-Davies predictions at late times and we explain why this should not come 
as a surprise. The mixed nature of the initial state in hyperbolic de Sitter space in fact 
ensures that all predictions match with those of the pure Bunch-Davies state in the planar 
de Sitter space that one started out with. The mixed nature of the initial hyperbolic vac¬ 
uum state therefore, in this case, does not imply observable departures from the standard 
planar Bunch-Davies predictions for the power- and bi-spectrum. Selecting instead the 
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pure hyperbolic vacuum does lead to differences that are however strongly suppressed in 
the curvature scale, as we again show explicitly by computing the power- and bi-spectrum 
for scalar field vacuum fluctuations. We end with a discussion on the implications in the 
context of (open) inflationary models and some remaining open questions that we hope to 
return to in future work. 

2 A family of de Sitter hyperbolic sections and their vacua 

We will start by constructing a limit in the family of de Sitter hyperbolic sections that 
reduces to the planar description of de Sitter space. This allows us to explicitly see the 
distinction between the natural vacuum state on a generic hyperbolic section and the 
Bunch-Davies vacuum state as defined using planar coordinates, which can be understood 
as a singular (but well-defined) limit of the natural hyperbolic vacuum. 

2.1 A generalized hyperbolic embedding 

Let us remind the reader that 4-dimensional de Sitter space can be dehned as the embedding 
surface in 5-dimensional flat space dehned by 

- + Xf + A| + A| + A| = 1 (2.1) 

where the de Sitter curvature length scale has been normalized to one. This embedding 
equation is clearly invariant under 50(1,4) transformations, corresponding to the isometry 
group of dS 4 . For our purposes we will be interested in two different coordinate sets 
on this embedding surface that both belong to the class of isotropic and homogeneous 
FRLW spaces. The standard coordinate set used for describing inhation is the planar 
one, identifying hat spatial sections. In what follows, we will suppress the two (spatial) 
coordinates X 3 and X 4 for purposes of efficiency, effectively suppressing an 5^ in the de 
Sitter space. The planar coordinates are dehned as 

Ao + = e**’ 

Xo - Ai = {rl - 1 ) ( 2 . 2 ) 

X 2 = Tp 

leading to the well-known planar expression for the induced de Sitter metric 

= —dtp + [drp + dn^] (2-3) 

with —00 < tp < - 1-00 and r > 0, and we reinserted the 5^ part. Since Aq -|- Ai > 0 these 
coordinates only cover the upper half diagonal part in the Aq versus Ai plane. Besides 
the obvious 50(3) isometries, the boost symmetries of the embedding space are realized 
on the planar metric as an isometry involving a particular combination of time translation 
and spatial scaling^. Because inhation redshifts away any existing spatial curvature present 
initially, this coordinate set should be an excellent approximation to derive the late-time 

^More precisely, it corresponds with the isometry t ^ t + 'y and r —>■ e 'r of the planar de Sitter metric. 
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effects of a sustained phase of cosmological inflation. Nevertheless, one could imagine 
a situation where our universe has originated from a tunneling event out of an eternally 
inflating false vacuum^. The nucleated bubble would have negatively curved spatial sections 
[16, 17], leading to the hyperbolic coordinate set 

Xq + Xi = cosh th + sinh th cosh rh 

Xq — Xi = — coshth + sinht/i coshr/^ (2.4) 

X 2 = sinh th sinh rh ■ 

With these identihcations the induced hyperbolic de Sitter metric reads 

= —dtf^ + sinht/i^ [dr^ + sinh^ rh dbl^] (2-5) 

where 0 < th < +00 and rh > 0 , and as before we included the full 5 ^ that was left out 
in the embedding coordinate identification. The coordinate singularity at th = 0 can be 
interpreted in the context of false vacuum decay as the creation of the open inflationary 
bubble. Note that th = 0 corresponds to Xq = 0 (and Xi = 1, X 2 = 0): the bubble 
nucleation time from the point of view of the embedding space. The spatial sections cor¬ 
respond to constant negative curvature slices that exhibit an S'0(l,3) isometry. We will 
in fact be interested in a one-parameter generalization of this hyperbolic coordinate em¬ 
bedding, obtained by boosting in the Xq-Xi plane of the 5-dimensional embedding space. 
Combined with rotations these transformation allow one to move the ‘nucleation’ time of 
the hyperbolic bubble to any specific point on the embedding surface. Just performing a 
Lorentz boost in the Xq-Xi plane will change the nucleation time (and position in Xi), 
which yields the following generalized hyperbolic coordinate set 

Xq -|- = e~'^ [cosht/j -|- sinht/j coshr/j] 

Xq — Xi = [—cosht/i -|- sinht/i coshr/i] (2.6) 

X 2 = sinht/j sinhr/i 

where 7 is the boost parameter. This generalized hyperbolic solution of the embedding 
equation will of course lead to the same induced metric, but the nucleation time and position 
of the associated bubble in the embedding space have now shifted to Xq = — sinh 7 and 
Xi = cosh 7 respectively. Moreover, since th > 0 one finds that Xq + Xi > e~'^, restricting 
the hyperbolic section to the upper right diagonal part in the Xq versus Xi plane, which 
overlaps with, but for any hnite 7 is smaller than, the part of de Sitter covered by planar 
coordinates. This is depicted in hgure 1. One can verify that in the limit of infinite 7 
the planar and hyperbolic coordinates cover the same region of de Sitter space, which is 
consistent with the observation that in this limit the hyperbolic nucleation time in the 
embedding space is shifted to Xq —)• — 00 . 

The generalization of the hyperbolic coordinate set introduced above allows us to 
explicitly relate the planar and hyperbolic sections of de Sitter space. Since the two co¬ 
ordinate sets cover the same region in the 7 — 7-00 limit, there should exist a one-to-one 

^Moreover, in the nineties models of open inflation were of particular interest, independent of whether 
their origin was due to tunneling [18, 19]. 
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Figure 1. Conformal diagram of dS 4 with the left- and right-hyperbolic patch as the upper-left 
resp. upper-right triangles. The dashed line is the unboosted situation 7 = 0 . For finite 7 (solid 
line), we see that that nucleation time of the left bubble gets pushed to earlier times, and vice versa 
for the right bubble. In the limit of 7 —>■ 00 , we can see that the left bubble will cover the entire 
upper-left triangle of the conformal diagram, coinciding with the planar patch. 


mapping between the coordinates in that limit. More precisely, we would like to introduce 
a new set of hyperbolic coordinates that are to be kept fixed in the limit 7 —)• 00 , and that 
in the limit exactly reproduce the planar coordinate embedding solution. Note that any 
(constant) shift or rescaling of the hyperbolic embedding coordinates is still a solution of 
the embedding equation, but will change the expression for the induced metric. Since we 
expect the range of the hyperbolic time coordinate to be extended to —00 and the negative 
curvature to be scaled away, we redefine 


4 = 4-7; 4 = 2 ^ he '^. 


(2.7) 


This leaves us with the following generalized hyperbolic solution to the embedding equation 


Xq + Xi = e [cosh {ih + 7 ) + sinh {ih + 7 ) cosh (24 e '’')] 

Xq — Xi = [cosh (4 + 7 ) — sinh (4 + 7 ) cosh ( 24 e”"’')] (2-8) 

X 2 = sinh {ih + 7 ) sinh (24 e~'^) 


where —7 < th < -|-oo. For finite 7 the shift in hyperbolic time and the rescaling of 
the hyperbolic radius (or equivalently the inverse rescaling of hyperbolic momentum) does 
obviously not affect any hyperbolic patch observables, but it does allow one to analyze the 
infinite boost limit in a simple and useful way. The induced hyperbolic metric now reads 


= —di^ -b sinh (4 + 7 ) 4e df\ -|- sinh (24e "*') dO. 


(2.9) 
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In the second line we performed the limit 7 —?• 00, keeping th and fh fixed, showing that ( 2 . 9 ) 
exactly reduces to the planar embedding solution ( 2 . 3 ). Note that all the 7 dependence in 
the induced metric is removed and one is left with precisely the planar line-element ( 2 . 3 ) 
in terms of the coordinates and r^. 

For any finite boost parameter 7 global de Sitter space is covered by two (adjacent) 
hyperbolic sections, see figure 1 . The other hyperbolic embedding can be obtained by 
changing the sign of Xi, resulting in the interchange of the expressions for Xq + Xi and 
Xq — Xi in ( 2 . 5 ). Acting with the same boost on this second hyperbolic embedding results 
in the opposite effect, moving the nucleation time to Xq —)■ +00. The opposite minus 
infinity boost should instead reduce to another planar section (with Xq + Xi and Xq — Xi 
in ( 2 . 3 ) interchanged), suggesting that the redefined coordinates in this case should read 

ih^th + i ; rh = ^rhe~'^. ( 2 . 10 ) 

Putting this together we obtain for the adjacent hyperbolic section the following generalized 
embedding 


Xq + Xi = —e [cosh {ih — 7 ) — sinh — 7 ) cosh { 2 fh e^)\ 

Xq- Xi = [cosh {ih - 7 ) + sinh (4 - 7) cosh { 2 fh e^)] ( 2 - 11 ) 

X2 = sinh {ih — 7 ) sinh (2f/j e^) 

where now <th < +00. The induced hyperbolic metric in this case is obtained by just 
replacing 7 with —7 in ( 2 . 9 ). By construction the limit 7 —>■ 00 should instead collapse 
and in a sense remove the adjacent hyperbolic section. Clearly, in the opposite 7 —)> —00 
the roles of the two hyperbolic sections are reversed. 

Having established this explicit relation between hyperbolic and planar coordinates, 
we can now use it to better understand and connect their respective vacua, which should be 
different for any finite value of 7. In particular, the planar Bunch-Davies state is known to 
be equivalent to the unique and de Sitter invariant Euclidean vacuum^. On the other hand, 
any pure hyperbolic vacuum state is defined on a negatively curved spatial slice that is not 
a de Sitter Cauchy surface. This means that the Bunch-Davies state in a single hyperbolic 
patch can only be described by an appropriately defined mixed state; see figures 2 and 
3 . The mixed state defined on one of the two (conjugate) hyperbolic sections reproducing 
the Bunch-Davies state was first constructed in [ 20 ] and was subsequently used in [ 21 ] 
to compute the reduced density matrix and the corresponding entanglement entropy for a 
single hyperbolic section. 

One application of the one-parameter family of hyperbolic de Sitter foliations is that 
one can confirm that the natural choice for a hyperbolic vacuum reduces to the planar 
Bunch-Davies state in the limit 7 —?• 00. Secondly, one could attempt to generalize the 
entangled expression for wavefunctions of the Bunch-Davies state, with support on both 
the left and right hyperbolic section, and work out its dependence on the embedding boost 

®The invariance of the Bunch-Davies vacuum under de Sitter isometries strictly speaking fails for mass¬ 
less fields, but since this subtlety does not affect our results we will ignore it from now on. 
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Figure 2. Conformal diagram of (iS '4 with the 
left- and right-hyperbolic patch shown. As nei¬ 
ther patch contains a Cauchy slice of the full 
dS4, restricting the Bunch-Davies vacuum to 
one of them will yield a mixed state. 


Figure 3. An observer confined to live in the 
left-hyperbolic patch (a bubble universe) can 
define his own pure hyperbolic vacuum. This 
state will differ significantly from the mixed 
state resulting from a restriction of the Bunch- 
Davies state to this bubble. 


parameter. In the 7 — 7-00 limit this should reduce to the pure planar Bunch-Davies state, 
implying that the reduced density matrix carries some non-trivial 7 dependence to make 
sure the associated entanglement entropy vanishes in the strict 7 —>■ 00 limit. 

To summarize, we established that the infinite boost limit of a hyperbolic de Sitter 
patch (and as a consequence also its corresponding vacuum state) reduces to the planar de 
Sitter patch (and the Bunch-Davies vacuum). This appears to be similar to an observation 
made in [ 22 ] where the static vacuum, understood as the empty state for a corresponding 
free-falling observer, was also argued to reduce to the Bunch-Davies state in the infinite 
boost limit. Note that to each hyperbolic patch one can associate a free-falling observer 
in one of the two center regions in between the hyperbolic patches that never intersects 
either one of them. These time-like curves are indeed connected to each other by the same 
embedding space boosts [23]. To complete the argument one needs to conhrm that the 
static vacuum state associated to this free-falling observer is connected to the hyperbolic 
vacuum state. Note that (for 7 = 0 ) the center region in between the hyperbolic patches is 
usually covered by coordinates that are obtained from the hyperbolic coordinates as follows 
th = iitc — f) and = re + if, resulting in the following center region metric 

ds^ = dtQ -b cos tc^ [~drc + coshre^dD^] , ( 2 . 12 ) 

where rc is now a time-like coordinate. Each of the two center regions clearly identifies 
a causal diamond belonging to the free-falling observer of interest. To make this explicit, 
one notices that the coordinate transformation r^ = sin tc and ts = rc indeed reproduces 
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the static patch metric, upon ignoring the two-dimensional sphere'^. This establishes a 
map from the hyperbolic to the static patch of a specific free-falling observer, relating their 
respective vacua and their behavior in the infinite (embedding) boost limit. 


2.2 The hyperbolic vacuum 

After having established a limit to obtain the planar embedding and coordinates, let us 
now remind the reader of the standard positive frequency modes on a single hyperbolic 
patch [20], as if it were the entire universe (see figure 3). The scalar wave equation for the 
hyperbolic patch of de Sitter (2.5) reads 


dt sinh^ tdt - 


sinh^ t 


sinh^ t 

where we defined the Laplacian on the three-hyperboloid 


= 0 


V^3 = 


1 


-9r.(sinh^ rdr) + — 


1 




sinh^ r sinh^r 


A natural set of solutions to the hyperbolic equations of motion (2.13) is given by 


sinh(t) 


(2.13) 


(2.14) 


(2.15) 


where it is customary to define u = 


I — The quantum numbers l,m label the usual 
SO{3) irreps, and together with the continuous quantum number p it completely specifies 
the hyperbolic momentum. Furthermore, i are the associated Legendre functions of the 
second kind and the are the orthonormal eigenfunctions of the hyperbolic Laplacian 
(2.14) 

V^3Tp;m(r,0) = -{l+p^)Ypim{r,VL). (2.16) 


For u > ^, there is in fact a supplementary set of solutions with p = — [20]. These so- 

called “supercurvature modes” will not be of interest for the purposes that are considered 
here, where our main focus will be on potential signatures in the large (“subcurvature”) 
momentum limit. We refer to [24] for an interesting account on the role and interpretation 
of these supercurvature modes. 

Switching to conformal time tj we can write the metric of the hyperbolic slice as 


ds^ = sinh^(t(r/)) {—drf' + dr"^ + sinh^(r)(in 2 ) ) (2-17) 

where p = ln(tanh(|)), or equivalently cosht = — ^ and —oo < r/ < 0. In terms of the 

conformal time p we find that in the far past p —)> —oo and in the limit of large momenta 
p ^ 1 one obtains 

pip /-1 \ ^ ^-ipp I'^ -^ A ^ ^-ipp^ ^2.18) 

\ tanhr/y \ ptanh py 

^That the part does not reproduce the standard static patch expression can be understood by realizing 
that this static patch region is rotated by an angle 7r/2 with respect to the standard embedding. This affects 
the angles, which have to be transformed as well in order to obtain the complete static patch metric. 
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so we can identify these mode functions with the “natural hyperbolic vacuum”: they define 
a state that is empty in the far past for large momenta, approaching the standard vacuum 
description in flat space. As expected, in the limit 7 —>• oo that we introduced in the 
previous section (2.10) the mode functions reduce to the standard Bunch-Davies mode 
functions in flat slicing, explicitly connecting the hyperbolic and planar patch vacua in this 
limit 

(2.19) 

For all the details we refer the reader to the appendix A.3, but it should be clear that 
the tildes on the coordinates in the above equation relate to the redefined hyperbolic 
coordinates that are kept fixed in the infinite boost limit. The mode functions (2.15) must 
of course be properly normalized, enforcing [bpim,bpim] — bii'6mm'b{p — p'), implying the 
following Klein-Gordon inner product 


lim P^^(cosh(f + 7 )) oc e I 1-^ 


7^00 


prj 


KG — bll' 6 mm'b{p P ) 


( 2 . 20 ) 


giving (see appendix A. 2) 


iV|,p = {P^p,P^p)kg 

2 sinh(7rp) 


( 2 . 21 ) 


TT 


With the help of (2.21) we can now express the field operator in a single hyperbolic 
patch as (keeping in mind that we are ignoring supercurvature modes) 


CO I 


/•OO ^ ^ 

JO r7 


1 1 
Npp sinh(t) 


bpimP"^_i (cosh(t)) Ypimir, n) + h.c. ) (2.22) 


/=0 m=—l 

defining the natural hyperbolic vacuum state as 

bpim I) — 0 '^Pi m. 


(2.23) 


This hyperbolic vacuum state can be understood as a natural choice in an isolated 
(stand-alone) open inflationary universe, as it is empty in the far past and reduces to the 
planar Bunch-Davies vacuum state in the infinite boost limit. As we have elaborated upon 
in the previous section, because a similar statement can be made for the vacuum in a static 
patch [22], and because the region in between the two hyperbolic patches contains a causal 
diamond region, this state can be analytically continued to the static (empty) vacuum for 
a free-falling observer that never intersects the two adjacent hyperbolic de Sitter patches. 

Clearly this state is very different from the unique de Sitter invariant Bunch-Davies 
vacuum for generic 7, so one would expect anomalous behavior similar to what happens 
in the de Sitter static vacuum or the flat Rindler vacuum. To that end let us analyze the 
behavior of the energy momentum tensor in a (generic) hyperbolic vacuum. Note that for 
a flat Rindler wedge in lightcone coordinates {u,v), there is a horizon at u = 0 and the 
Fulling-Rindler vacuum \0fr) corresponds to the empty state in a single wedge. In that 
case it is well-known that the T^u component of the energy momentum tensor (with the 
usual UV-divergence removed by subtracting the UV-divergent expectation value of in 
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the Minkowski vacuum |0 m)) diverges as one approaches the horizon: {Tuu)fr — {Tuu)m = 
— in 1 + 1 dimensions for u > 0 (for a nice derivation of this result see [25] ). 

A similar analysis can be done for the energy momentum tensor in the hyperbolic de 
Sitter patch, where the global de Sitter invariant vacuum state is now the Bunch-Davies 
vacuum The obvious difference with the Rindler wedge is the absence of a timelike 

Killing vector. In addition, the t = 0 surface is a (light-) cone, so a better analogy is 
with Milne space, to which the de Sitter hyperbolic section reduces for small t. In any 
case, we will use the same regularization procedure, restricting to the minimally coupled 
massless case i/ = |. The most convenient method to calculate components of the energy 
momentum tensor makes use of the Wightman function G~^{x,x 'and specifically we 
will look at the following contribution 

{{da4>)‘^{x,t)) = lim dadaG{x,x',t,t'). (2.24) 

The Wightman function for the Bunch-Davies state is well known, but here we use the 
expression in terms of an integral over the hyperbolic momentum p as given in [20] . This 
allows us to consistently regulate the UV-divergence of in the two states of interest. 
In appendix C we show that the difference {Ttt)H — {Ttt)BD is UV-finite and diverges as 
t 0 

{Tu)h - {Tu)bd = + ^ (^) • 

So we conclude that the hyperbolic vacuum JD//) has singular properties that are completely 
analogous to the Minkowski Fulling-Rindler, Milne and de Sitter static vacuum, see also 
[26]. The energy momentum tensor diverges in the limit t —)• 0, so infinite energy seems 
to be required to prepare the state at the (singular) origin. A complete description all the 
way until t = 0 is therefore obviously inconsistent, but strictly speaking that does not need 
to be fatal in a cosmological setting, in the sense that in a stand-alone open universe this 
might be interpreted as the Big Bang singularity. 

Of course, arguably the most natural and well-behaved choice for an initial state on 
hyperbolic de Sitter sections is the de Sitter invariant Bunch-Davies state, to which we 
turn next. 


2.3 The Bunch-Davies state in the hyperbolic patch 

Here we will just briefly summarize the results of [20] and [21]. More details can be found 
in those papers and in appendix A. The most important observation is that mode functions 
of one of the hyperbolic patches (2.15) do not correspond to regular mode functions on 
the full (Euclidean) de Sitter space. In [20] the hyperbolic mode functions are analytically 
continued to the other hyperbolic patch, allowing them to construct a set of regular mode 
functions that can cover all of de Sitter space as follows 




isin(7r(i/-|)) ip , . ism{-n{ip+v-\))eT[v+\+ip\ ip , s 

sinh(p7r) sinh(7rp) Y[v+\-ipy ' 


for z G i? 
for z G L 


(2.26) 
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These mode functions do not yet describe the Euclidean or Bunch-Davies vacuum, which 
can for instance be concluded by the fact that they are not (anti-)symmetric under the 
transformation R ^ L. It turns out that the linear combinations XR XL correspond 
to the proper mode functions associated with the Euclidean or Bunch-Davies vacuum, as 
was proven by computing the Wightman function [20]. The (still to be normalized) mode 
functions are linear combinations of the associated Legendre functions 


Xp,cr — 



(2.27) 


where a = ±1^ and z = cosh(t); the expressions for the a’s and /3’s are given in (A.7). We 
stress that the associated Legendre functions P in (2.27) do not have to be analytically 
continued any further®. The full field expansion, with creation and annihilation operators 
a„pini satisfying [dapZm, 0^/^///^'] = - p'), IS given by: 


(t){t,r,n) 



{^aplmXp,o-i,^)^plm{L, 12) T h.C.) 


(2.28) 


where is the Klein-Gordon norm consistent with the commutation relations'. In 

appendix A.2 we show that the normalization Ny.pt7 is given by 

R^^pty — {Xp,cr^plmj Xp,cr^plm,)KG 


q=L,R 


(2.29) 


where a, f3 denote the complex conjugates of a, (5. We conclude that the Bunch-Davies 
vacuum state, in a single hyperbolic patch, is defined as® 


aapim\^BD) = ^ Vcr,p,/,m. (2.30) 

Now let us describe the relation between the creation and annihilation operators of the 
modes (2.27) and the creation and annihilation operators of hyperbolic modes (2.15). Both 
field expansions (2.28) and (2.22) are linear combinations of associated Legendre functions. 
We can hnd the relation between the do-pim and the bqpim {q = L, R) by comparing the 
coefficients 

E IXpp ( ^ ~ X 

~jX \^p,qb'aplm + l^p,q^apl-m ) ' (2-31) 

a=±l 

Given (2.29,2.31), they enforce 


Scrp/m, 

— 5{j) p 


✓ 

bqplm^ bg'pfPrn' 

5(^p p 5iif dqq 



> 4^ < 





= 0 



bqplrri'; bqfp'l'm' 

= 0 


(2.32) 

= ±1 is related to the combination ± 

®They are constituents of xl and Xfl, which are already regular everywhere. This is different from [13]. 
^Strictly speaking the expression (2.28) is incomplete, since we should also include the “zero mode”. For 
our purposes this will however not affect the results. 

®As before we ignore the supercurvature modes. 
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For more details we refer to [20] and the appendices. The above relationship con¬ 
firms that the normalizations (2.21) and (2.29) are consistent and in particular that the 
right normalization for the hyperbolic mode functions is given by (2.21). This will be of 
importance when comparing the predictions for the power- and bi-spectrum of the two 
different states under consideration: the pure hyperbolic vacuum and the Bunch-Davies 
state (as we will do in section 3). The latter is a mixed state from the point of view of a 
single hyperbolic patch, due to the entanglement between the two hyperbolic patches in 
the Bunch-Davies vacuum. 

Let us here remind the reader that we would like to compare the predictions for the 
expectation values of (scalar field) quantum fluctuations in the two different states that were 
introduced above. A priori different initial states give different predictions for the cosmic 
microwave background temperature anisotropies and the large scale structure distribution. 
We should stress that we are technically not considering an actual bubble nucleation event, 
where more intricate and model-dependent bubble wall physics could lead to additional 
effects [27, 28], see also [29]. Instead, we will work under the assumption that the two 
states that were introduced capture an essential difference that is generic: the entangled 
nature of the Bunch-Davies vacuum implies a mixed initial state, whereas the hyperbolic 
vacuum corresponds to a pure state on a single hyperbolic section. Different (presumably 
more realistic) states in open universes have been considered in the past [14, 27], but these 
states appear to be pure (and excited) hyperbolic states, so they do not capture any effects 
related to the mixed nature of the initial state. In the process we hope to clear up some 
confusion that might have arisen and that could also have consequences for more realistic 
bubble states that were considered in the past. 

We should add that one might anticipate the differences between the two states to 
only become visible at small hyperbolic momentum p < 1, i.e. scales comparable to the 
hyperbolic curvature. However, even small curvature suppressed changes in the initial state 
might be enhanced in the (nonlinear) bi-spectrum, as has been pointed out and analyzed in 
[30, 31] and for a certain generic type of mixed state in [15]. This motivates our particular 
interest in computing the bi-spectrum, comparing the different states to the planar Bunch- 
Davies result. But hrst let us review some general facts regarding correlators in de Sitter 
space and summarize the results for the two-point functions. 

3 Correlators in hyperbolic de Sitter space 

Making use of the previously established relations between the de Sitter invariant Bunch- 
Davies vacuum and the hyperbolic vacuum state, we will compute both the Bunch-Davies 
and the hyperbolic vacuum power-spectra of scalar held quantum huctuations®. The 
Bunch-Davies result can also be calculated using a reduced density matrix formalism in the 
hyperbolic patch. Let us from the outset emphasize that within our basic de Sitter set-up, 
even though the Bunch-Davies state is mixed from the hyperbolic patch perspective, all hy¬ 
perbolic Bunch-Davies correlators should match the (hyperbolic coordinate transformed) 
planar Bunch-Davies correlators. As a consequence one can rule out large deviations of 

®See [26] for related work on the response of Unruh detectors. 
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Bunch-Davies hyperbolic correlators at late time and large momenta (when the hyperbolic 
coordinates reduce to planar coordinates) as compared to the planar Bunch-Davies corre¬ 
lators. That leaves the (pure) hyperbolic vacuum as the potentially more interesting state 
to consider, as far as enhanced initial state effects with respect to the planar Bunch-Davies 
state are concerned. 

Let us start by pointing out that the field operator (pp evaluated on points in the left 
hyperbolic patch is trivial in the right hyperbolic patch 


(f)p{x) = (j)L,pix) (8) Ir for X G L. 


(3.1) 


As a consequence, de Sitter n-point functions of fields (/)p in the Bunch-Davies state, eval¬ 
uated on points in the left hyperbolic patch, can be calculated either using the full global 
description or by using a reduced density matrix pL = Tr-^^ {|DBD)(flosol}’ their 
results should agree. This is shown explicitly in appendix B.2. By defining the bqpim as in 
(2.31), we can write the field operator for arbitrary values of p, I, m as 


4^plm{x) — b^pijn 


1 


+ bRplm r^plm + h-c. 


'Np: 


(3.2) 


where 



1 (t) for t G L 
0 for t G P 


(3.3) 


and vice versa for Although these functions are not mode functions on a full 

U— 2 

Cauchy slice covering the de Sitter space, we are allowed to express the field in terms of 
them. Note that the bi and bp operators mutually commute. To make explicit that the 
field operator decomposes in the left and right hyperbolic patches, we write 


*/’pZm(^) — I bpplm 


J—PT-1 T^plm + h-C. ) (8) liJ 

rspp ^ 2'^ 


1 


+ II ® ( bRplm Pf_l^p^plm + h-C. 


(3.4) 


Note that the expansion of a scalar field in Minkowski spacetime in terms of left and right 
Rindler wedge modes is similar (see for instance [25]). The restriction of the operator (3.4) 
to points in the left hyperbolic patch is by definition equal to the full operator evaluated 
on points in the left hyperbolic patch. Note that if the field operator evaluated on points in 
the left hyperbolic patch would also have support on the right hyperbolic patch, it would 
not make sense to do a density matrix calculation as done above. 

Clearly therefore Bunch-Davies scalar field correlators should be the same, indepen¬ 
dent of whether one uses hyperbolic or planar coordinates. Of course, since the de Sitter 
invariant length is expressed differently in terms of planar or hyperbolic coordinates, the 
functional dependence of the equal (hyperbolic) time correlators will look different. Since 
the difference between planar and hyperbolic coordinates vanishes in the late time and 
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large momentum limit, the the planar and the hyperbolic Bunch-Davies correlators match 
in that limit and small modifications are suppressed in the hyperbolic curvature scale. We 
conclude that hyperbolic Bunch-Davies correlators can be computed either using a global 
de Sitter description (for which the Bunch-Davies state is a pure initial state) or by con¬ 
sidering a single de Sitter hyperbolic patch (for which the Bunch-Davies state equals a 
mixed initial state). 

After these important preliminaries let us now proceed by computing the power spec¬ 
trum of a massless scalar field in a hyperbolic coordinate patch in the hyperbolic vacuum 
and Bunch-Davies initial state respectively, as a function of the hyperbolic momentum p. 


3.1 Power—spectrum results 


For most of the details we refer to the appendix B. Here we will just quote the main results. 
For the two point function in the hyperbolic vacuum state we find 


{Qh\4>p(I>'p\^h) = S{p-p'] 


p cosh^ (t) + p^ 
sinh^(t) 47r2 p^ + I 


(3.5) 


where we have used the completeness relation of the eigenfunctions of the hyperbolic Lapla- 
cian^*^ and the commutation relations. At late times t —)> oo this approaches 




H‘^p 


P P -|- 1) 

and the appropriately normalized power spectrum (at late times) equals 


kH 


ip) = 




p- 




(3.6) 


(3.7) 


The same (hyperbolic coordinate patch) two point function in the Bunch-Davies vac¬ 
uum is instead found to be equal to 


{0,bd\4>p4’p'\^bd) = ^{p — p) — 




p cosh^ (t)-|- 


sinh^ (t) dvr^ p^ + 1 


coth( 7 rp). 


(3.8) 


This result can either be obtained from a direct calculation using the global Bunch-Davies 
vacuum construction and restricting to one of the hyperbolic coordinate patches [20], or 
from a (mixed) density matrix calculation in a single hyperbolic coordinate patch, using 
the explicit expression for the density matrix as reported in [21], as we confirm in appendix 
B.2. As alluded to earlier, the reason for this expression to not exactly reproduce the scale- 
invariant planar coordinate result for the scalar field power spectrum in the Bunch-Davies 
vacuum is that different coordinates are used. As the hyperbolic and planar coordinates 
are the same at late times and for small distances, the late-time power spectra at large 
momentum should be the same as well. At late times t —)• oo we find 


i^BD\(l>p(l>p'\^BD) 


H^p 


d7r2(p2 _j_ 


coth( 7 rp). 


10 


|lplm| — 2^2 • 


(3.9) 
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(3.10) 


Correspondingly, the power spectrum (at late times) is given by 


A 


,BD 


(P) = 




P 




coth(7rp). 


A"(P) 



Figure 4. The power spectra (logarithmic scale) for the hyperbolic vacuum (blue) and the Bunch- 
Davies vacuum (red), as function of the hyperbolic momentum p with H = 1. The dashed line 

tt2 

indicates the scale-invariant planar Bunch-Davies result 

Looking at these power spectra we indeed find that for p ^ 1, when coth vrp ~ 1 and 
+ 1 ~ , the hyperbolic vacuum as well as the hyperbolic Bunch-Davies result matches 

the standard scale invariant planar Bunch-Davies result They only start to differ 

from each other and the standard planar Bunch-Davies expression for sufficiently small 
momenta p ^ 1 (see figure (4)). Note that although the corresponding wavelengths are 
expected to lie far outside our observable window, given the fact that they correspond to 
length scales longer or comparable to the hyperbolic curvature scale, for both hyperbolic 
states the power is suppressed as compared to the standard planar Bunch-Davies result. 
As these departures from the standard planar result become evident, one should keep in 
mind that the difference found in the hyperbolic Bunch-Davies result can be attributed 
to a coordinate change, whereas (part of) the change in the hyperbolic vacuum power 
spectrum is related to an initial state modification. Although this difference might seem 
unimportant at this point, when considering the bi-spectrum in the next section it is a 
relevant distinction, since the bi-spectrum has been found to be particularly sensitive to 
changes in the initial state. 

Finally let us also remark that the bubble state put forward in [27] and more recently 
used in [14] does resemble the Bunch-Davies state in the sense that the expectation value 
of the number operator on hyperbolic sections agrees (giving thermal occupation numbers 
in terms of co-moving momentum), but in one aspect is crucially different due to the fact 
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that it has been constructed as a Bogoliubov transformation of the hyperbolic vacuum. As 
a consequence the state necessarily corresponds to a pure state (and inherits the singular 
properties of the hyperbolic vacuum) and it explains why their result for the power spec¬ 
trum does not agree with our hyperbolic Bunch-Davies result that necessarily corresponds 
to a mixed state. 

We conclude that, independent of the particular initial state under consideration, any 
power-spectrum signatures of an open inflationary universe are confined to the curvature 
scale, which has to be several orders of magnitude larger than the largest observable length 
scale in the universe. Although the initial hyperbolic state is mixed when assuming a 
(globally defined) planar Bunch-Davies state, the power-spectrum results in this admit¬ 
tedly basic set-up in which all bubble wall physics is ignored, do not show large deviations, 
as should be expected. In fact, this statement is true for general hyperbolic n-point correla¬ 
tors in the Bunch-Davies vacuum. Potentially enhanced bi-spectrum results due to initial 
state excitations, as compared to the standard planar Bunch-Davies result, can therefore 
only be expected assuming the pure hyperbolic vacuum as the initial state in the hyperbolic 
patch. 

3.2 The bi—spectrum in the hyperbolic vacuum 

Let us next consider the bi-spectrum of scalar density perturbations in the hyperbolic 
vacuum, which can be viewed (for any finite boost parameter 7) as some sort of excited 
initial state with respect to the standard planar Bunch-Davies vacuum. We will mainly 
be interested in the so-called squeezed limit, for which previous work uncovered enhanced 
results for excited (planar) Bunch-Davies intial states [15, 30-32]. Moreover, as before it 
should be reasonable to work in the sub-curvature approximation p 3> 1, for which the 
hyperbolic momenta pi are approximately equal to the standard (flat) wavenumbers ki up 
to curvature suppressed corrections. 

To compute the hyperbolic bi-spectrum we need the action to third order in the scalar 
density perturbation cp and the hyperbolic curvature introduces some new ingredients as 
compared to the planar calculation [33], which have been carefully dealt with in [14]. The 
state of interest in [14] is however the bubble state constructed first in [27], which is related 
to the hyperbolic vacuum by means of a Bogoliubov transformation and crucially differs 
from the Bunch-Davies vacuum. The bubble state can be viewed as a (well-motivated 
and first principles derived) initial state modification with respect to the hyperbolic vac¬ 
uum, explaining their interest in trying to identify enhanced features in the bi-spectrum. 
Note that their computation is not applicable to initial state modifications of the (planar) 
Bunch-Davies vacuum. Such an interpretation would only be valid in the infinite boost 
limit, when the hyperbolic vacuum reduces to the planar Bunch-Davies vacuum. As we 
explained, in our pure de Sitter set-up it is the hyperbolic vacuum initial state that (for 
finite boost parameter 7) corresponds to an excited state as compared to the planar Bunch- 
Davies vacuum and could therefore potentially display interesting bi-spectrum enhance¬ 
ments. Instead the hyperbolic Bunch-Davies bi-spectrum equals the planar Bunch-Davies 
bi-spectrum and all apparent changes can be related to the coordinate change from planar 
to hyperbolic. So here we will be interested in computing the results for bi-spectrum in the 
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hyperbolic vacuum, which can fortunately be extracted straightforwardly from the results 
in [14], 

So let us first briefly review the basic results reported in [14] and then apply them to 
our case of interest. We will be interested in a massless minimally coupled scalar field, for 
which the positive and negative frequency modes dehned in a single hyperbolic patch, can 
be nicely expressed as 


^ ^cosh^-ipsmhnj„ 

y2p(l+p2) 

The perturbed metric in one of the hyperbolic patches, is written in the ADM formalism 


as 


ds'^ = -N^dt^ + hij {dx^ + N^dt) {dx^ + NUt) , (3.12) 

where, as usual, N is the lapse function, W is the shift and = a?{t)e^^^ij is the spatial 
metric, with curvature perturbation Q. We will be interested in gravitationally induced non- 
linearities on the scalar density perturbation, requiring that we need to introduce (generic) 
slow-roll evolution of the background scalar (^{t) in order to couple the scalar inflaton field 
to the scalar density perturbation. Now, plugging the above metric into the action for the 
scalar degree of freedom (assuming slow-roll evolution) and solving the constraint equations 
order by order, one can obtain the quadratic and cubic (and higher) action. The quadratic 
action for scalar perturbations, in the flat gauge C = 0, to leading order in the slow roll 
parameters is 

<5^^^ = j dtd^xa{tf^ , (3-13) 

where di is the covariant derivative with respect to Up to this point this should all be 
familiar, so let us now turn to the cubic action on a hyperbolic patch. As explained in [14] 
the dominant term in the third order Lagrangian is found to be 

£(3) = (^(a^ - 3)"^ (3-14) 

The new ingredient due to the hyperbolic curvature in this action is the —3 term. The 
is (as usual) the redefined field (^ —)• £c), dehned as 


£ = £c + 


4(a/a) 


(5^ - 3) ^ dupcd'ipc , 


(3.15) 


which does not affect the quadratic action and has removed terms in the cubic action that 
are proportional to the equations of motion. The bi-spectra of £ and £c are then related 
as follows 


(£( xi )£( x 2 )£( x 3 )) = {(Pc{xi)ipc{x2)^c{X‘i)) 


+ 


4(d/a) 


(9^ - 3) {diipc{xi)(pc{x 2 )) {d\c{xi)^Pc{x'i)) + permutations 

(3.16) 
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In harmonic space, introducing the geometrical factor 




(3.17) 


the bi-spectrum B{pi,p 2 ,P 3 ) is dehned as 

{'fpihm,i^P2l2m2fp3l3m3)=B{pi,p2,P3)Yp\p^p^Gl^l^l^^ (3.18) 


Of the two bi-spectrum contributions due to the field redefinition, the first term is then 
computed using eq. (3.14) in the in-in formalism, where we already assumed the sub¬ 
curvature limit allowing one to replace the momentum p with the planar momentum k 


B{ki, k 2 , fes) = 2Re 


ivki{0)vk2{0)vk3{0) 



^d)uk2 {v)uk3 (rj) 


-hlo2-hle-^3 . 

(3.19) 


In a hyperbolic (quasi) de Sitter space the integral over conformal time is naturally divided 
into two eras: from —oo to —1, where the curvature term is dominant and from — 1 to 0 
corresponding to standard inflationary expansion. Assuming slow-roll evolution in both 
eras, the expressions for (j) in terms of the slow-roll parameter e change when transitioning 
from the curvature-dominated era into the inflationary era. The mode functions defining 
the hyperbolic vacuum in the sub-curvature approximation p ^ 1, in the different 

eras, read as follows 


Ukin) - - 
Ukiv) - - 
Ukiv) = 


a(ry)\/^ 
H 

duk 1 


i-oo<v<-i) 


{l + ikp)e (-1 ^ < 0) 


1 


Vk{0) ^ 


dp a{p) a?{p) 

H 

vw 


-e 


—ikr} 


(3.20) 


Using these expressions to do the integral in (3.19)and ignoring effects due to the discon¬ 
tinuity in the transition between eras, one obtains [14] 


B{ki,k2,k3) = 




Akik2k3 [ki + k2 + k3) ykl k^ kU 


1 1 1 

+ To + 


(3.21) 


This result agrees with the standard planar Bunch-Davies result in single held inhation. 
For completion the additional contribution from the held redehnition equals 


B^^'^^^iki,k2,k3) = +Pe™^tations, (3.22) 

but as shown in [14] the contribution from this term is subdominant in the squeezed limit 
fes <C fei ~ A:2 = A: and therefore can be ignored when looking for other enhanced contribu¬ 
tions in the squeezed limit, as expected for initial state modihcations. The enhancement of 
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initial state modifications in the squeezed limit, as compared to the hyperbolic vacuum, was 
indeed observed after adding the specific negative frequency term to the mode functions 
that describes the bubble initial state [14]. However the Bogoliubov coefficient describing 
this excited bubble state is exponentially suppressed in the hyperbolic momentum over 
the curvature scale (which is normalized to one) e~'"^, meaning that these effects for the 
observable modes in the sky are exponentially suppressed and undetectable. 

After having reviewed the hyperbolic bi-spectrum calculation in [14], we will now 
connect these results to our case of interest. To start out with, the quoted result for the 
bi-spectrum (3.21) is of course the hyperbolic vacuum bi-spectrum in the late time, sub¬ 
curvature limit, which agrees with the standard planar Bunch-Davies result. Clearly then, 
no (modified initial state) enhancement with respect to the planar Bunch-Davies result is 
found in the late time and large momentum limit. The squeezed enhancement that was 
revealed in [14] for the excited bubble state, albeit exponentially suppressed in momentum, 
is with respect to the hyperbolic vacuum and can only be interpreted as enhancement with 
respect to the planar Bunch-Davies vacuum in the late time and large momentum limit. 
As we noted, the late time and large momentum limit effectively corresponds to the action 
of the inhnite boost in the embedding space for which the hyperbolic vacuum indeed 
matches the planar Bunch-Davies vacuum state, explaining the result. The excited bubble 
vacuum, obtained through a Bogoliubov rotation from the hyperbolic vacuum and therefore 
corresponding to a pure state, will also not agree with the hyperbolic Bunch-Davies result. 
Instead the hyperbolic bi-spectrum in the Bunch-Davies state has to match the planar 
Bunch-Davies result (written in hyperbolic coordinates), obviously excluding a (squeezed) 
enhancement. 

Our first conclusion is therefore that the bi-spectrum in the hyperbolic vacuum is 
clearly (and unsurprisingly) not enhanced with respect to the planar Bunch-Davies result 
in the late time and sub-curvature limit. In fact it is straightforward to (perturbatively) 
extend this conclusion beyond the strict sub-curvature or large momentum limit by ana¬ 
lyzing the leading correction. The first correction in the large momentum approximation 
is obtained from an asymptotic expansion of the hyperbolic mode function at late times 

Upiv <Cl,p > 1) « -^= (I -h ipn) ^ • (3-23) 

Noting the last term (in brackets) one observes that the leading large momentum correction 
will give additional contributions suppressed in the large momentum limit with at least one 
factor of but leaving the relative momentum dependence in tact. As a consequence 

enhancements in the squeezed momentum limit as compared to the planar Bunch-Davies 
result are excluded, as that would typically require additional terms featuring negative 
frequencies. It should be clear that corrections suppressed as 1/p^, without additional 
enhancements in particular momentum configurations, does not constitute an interesting 
non-Gaussian signature of an open inflationary universe. 

^^Re-installing the hyperbolic curvature scale, which was set to one, this term would be explicitly dimen¬ 
sionless and read fec/p^. 
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We conclude that, unfortunately, the bi-spectrum in the hyperbolic vacuum does not 
produce interesting enhancements that could be searched for. This conclusion is in fact 
corroborated by effectively constructing the hyperbolic vacuum as an excited state on top 
of the (global) planar Bunch-Davies vacuum (see appendix D). Effectively one then discov¬ 
ers that the relevant Bogoliubov coefficients are suppressed exponentially in momentum, 
ensuring that these effects will not be observable. As emphasized, the same is true for the 
bi-spectrum in the hyperbolic Bunch-Davies initial state. So even though the hyperbolic 
vacuum can be thought of as an excited state with respect to standard planar Bunch- 
Davies vacuum, it is not of a type that gives rise to large (enhanced) corrections in the 
bi-spectrum as compared to the standard planar Bunch-Davies result. 

4 Conclusions 

Before summarizing our results, let us remind the reader once more that our motivation was 
to carefully study the relation between the hyperbolic and planar coordinate patches and 
their corresponding states in (mostly) pure de Sitter space. We believe these results to be 
of interest, and partially applicable, in the context of de Sitter false vacuum decay, but it is 
also clear that in that case a more complete analysis should include (model-dependent) wall 
physics that will affect the details. Instead we concentrated on a general and qualitative 
difference between two examples of initial states on a hyperbolic section of de Sitter space: 
the pure hyperbolic vacuum and the (entangled) de Sitter invariant Bunch-Davies state, 
corresponding to a mixed state. We first of all noted that the pure hyperbolic vacuum is 
formally inconsistent, due to the energy momentum tensor becoming singular at the null 
boundary of the hyperbolic section. This issue should plague all pure hyperbolic states, in 
particular also the bubble state introduced in [27], which is obtained from the hyperbolic 
vacuum by a (unitary) Bogoliubov transformation. 

The Bunch-Davies vacuum is qualitatively different, giving rise to a mixed state on a 
single hyperbolic section. The fact that the Bunch-Davies state is mixed on a hyperbolic 
section could a priori imply interesting non-Gaussian signatures [13] as compared to the 
bi-spectrum in the hyperbolic vacuum. But as we pointed out all planar n-point correla¬ 
tors in the Bunch-Davies state, with all points in the hyperbolic section, should equal the 
hyperbolic n-point Bunch-Davies correlators. Any differences are generated by the coor¬ 
dinate change from planar to hyperbolic, instead of arising due to changes in the initial 
state. So even though the Bunch-Davies state is a mixed state on a hyperbolic section, 
certainly implying different results as compared to the hyperbolic vacuum, it will reproduce 
all planar Bunch-Davies correlators, up to coordinate changes. Since the hyperbolic and 
planar coordinates coincide in the late time and large momentum (sub-curvature) limit all 
effects due to the coordinate change should disappear and the leading corrections will not 
be enhanced. 

We stressed that instead the hyperbolic vacuum should generically be considered as 
an excited state as compared to the Bunch-Davies vacuum. As a consequence that should 
be the state of interest for computing the bi-spectrum and look for enhanced features. 
Using the hyperbolic coordinate embedding we explicitly constructed a family of hyperbolic 
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solutions that reduces to the planar coordinates in the inhnite boost limit, as such providing 
a limiting relation between the hyperbolic vacuum and the planar Bunch-Davies vacuum. 
As a corollary we also argued that the hyperbolic vacuum can be mapped to a specific static 
vacuum, implying that the static vacuum should also reduce to (a sector of) the Bunch- 
Davies state in the inhnite boost limit, as was hrst noted in [22]. Again, this limiting 
behavior implies that in the late time and large momentum limit, the bi-spectrum results 
for the hyperbolic vacuum should agree with the standard planar Bunch-Davies result and 
using the results of [14] this was indeed conhrmed. Looking at the leading correction in 
the large momentum expansion, we verihed that no enhancement in particular momentum 
conhgurations is generated and that the corrections are at least suppressed as 1/p^. These 
type of curvature suppressed Non-Gaussian corrections will clearly be impossible to detect. 
So unfortunately, on the basis of our analysis here, we should conclude that no detectable 
signals of an open inhationary universe in the huctuation statistics is expected on small 
sub-curvature scales. 

To summarize the two hyperbolic states introduced, one of them mixed (Bunch-Davies) 
and the other one pure (hyperbolic vacuum), make almost identical predictions in the late 
time sub-curvature limit. In fact, in the inhnite boost limit the states become formally iden¬ 
tical to the planar Bunch-Davies vacuum. For the mixed Bunch-Davies state this seems to 
imply that the density matrix pbd should depend on the boost parameter 7. Correspond¬ 
ingly, the associated von Neumann entropy Tr (—pbd Ih/Obd) of the mixed Bunch-Davies 
state on the hyperbolic section should depend on the boost parameter 7 to ensure that 
the entropy vanishes in the inhnite boost limit. This density matrix was computed in [21] 
and it would be of interest to consider the generalization for non-zero boost parameter 7. 
Although one might think the density matrix and corresponding entropy to be boost invari¬ 
ant, this is not entirely obvious and the above observation does indeed suggest it might not 
be, perhaps in some subtle (singular) way. The dependence on the boost parameter shonld 
be such that it is invariant under 7 —)> —7, effectively interchanging the two hyperbolic 
sections. Since the boost dependence can be implemented through a simple rescaling on 
the left hyperbolic momenta (and a time shift) 2.7 and the inverse rescaling on the right 
hyperbolic momenta (and time shift) it should be possible to trace the boost dependence of 
the Bunch-Davies state in terms of the left and right hyperbolic modes. It should then be 
straightforward to construct the corresponding density matrix and conhrm that the density 
matrix and corresponding entropy become trivial in the inhnite boost limit. 

In more realistic scenarios, trying to incorporate the bubble nucleation dynamics, an 
initial state has been proposed that can be constructed by performing a unitary Bogoliubov 
transformation on the hyperbolic vacuum state [27]. As a consequence this state appears 
to be pure, which at hrst sight contradicts the general expectation that the initial state 
inside a bubble, after false vacuum decay, should be mixed. Although the pure bubble 
state is in certain aspects similar to the Bunch-Davies vacuum, the correlation functions 
will be different as compared to the Bunch-Davies hyperbolic correlation functions. In 
addition, as in the hyperbolic vacuum, the expectation value of the stress tensor in this 
pure hyperbolic bubble state should become singular as one approaches the bubble wall. 
It might be of interest to revisit the original construction and see how it can be adjusted 
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to construct a mixed hyperbolic initial state instead. We hope to address this and some of 
the other remaining questions in future work. 
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A Mode fnnctions 

A.l Solutions to the hyperbolic equation of motion 

The metric for both the left and right hyperbolic patch is given by: 

^ (—+ sinh^ t [dr‘^ + sinh^ r dDl)) (A-l) 

where the coordinates are dimensionless and c = 1. The action for a massive 

non-interacting minimally coupled scalar field (p is given by: 

S = -^ j {g^''d^(j)dy4> + rn^cj)^) . (A.2) 

The action of a conformally coupled scalar field can be written in the Einstein frame with 
effective mass m? = 2H^. The equation of motion is given by 


0 = 


1 d n d 1 

-5—— smh (tj— -^— 

sinh^(f) dt dt sinh^(f) 


v: 


m 

Jp 


1 


1 


d . , 3. . d 

-o —smh it)— -o— 

sinh^(f) dt dt sinh^(f) 




(A.3) 


where u is defined as z/ = y | ^ and V^3 is the Laplacian on the hyperboloid T-L^. We 

will use u' = u — consistent with [20], such that u' = 1 corresponds to the massless 
minimally coupled case and = 0 corresponds to the massless conformally coupled case 
for which the effective mass is = 2H^. 

The eigenfunctions Ypi^ of the Laplacian on the hyperboloid that are regular 
in r = 0, are given by [20]: 


Y .^3 Ypljyi 

Yplmir 1 

fpl{r) 


^')Yplm 

)Ylmm 

+ ^ + _ P _p-'-1/2 . ^ p 

/lAhzp±tLsmh' ''' 


/ sinpr 


__ smh V _ I _ 

7 rT{—ip + l + l) d(coshr)* \sinhr 


(A.4) 
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where Yim{^) is the normalized spherical harmonic function on the unit two-sphere, r( 2 ;) 
is the Gamma function and -P^(z) is the associated Legendre function of the first kind [34]. 
The mode functions that correspond to the natural hyperbolic vacuum are given by; 

^ (cosht) positive energy modes ' 


sinht ‘ p ^ 0 

(cosht) negative energy modes ~ 


(A.5) 


Mode functions on a Cauchy slice of de Sitter, must be regular and consist of linear 
combinations of the hyperbolic mode functions in the left and right hyperbolic patches [20] . 
The mode functions that correspond to the Bunch-Davies state are given in [20]: 


Xp,t 


= 


^ p^P(^) _ I_ p- 

T{v' + ip+l) T{v' — ip + l) ''' 

rrp'^P — p-™'^' . n-p~'^P — p-™'^' 

K^i^)- Z > 

r(i/' — +1) 


-ip/ 


-ip/ 


T{u' + ip+l) 


for X G i? 


for X G L 


(A.6) 


^ I for X G ii 

1 + Pp,lP7^{^) for a; G L 


where z = cosht and the constants ^ and ^ are defined as: 


P:<J 


^p,R T(y' -\-ip+l) 


Pp,L = 




r{i''-ip+i) f^p,R 


^Ir = - 


r{i/'—ip-\-l) 


(A.7) 


These mode functions must be normalized through the Klein-Gordon normalization (see 
section A. 2). 

A.2 Klein—Gordon normalization 
Hyperbolic modes 

We normalize the hyperbolic modes on the hyperbolic patch using the variable z = cosht 
and using the orthonormality of the Ypim- 

^pp — {4^plm: 4^plm}KG 

= Z / dTi^ (^(l)pijndf^(f> — (f) 

Jy. 

= i sinh^ t { _ -P;;'^(cosht) ^^ / Pj(cosht) \ \ (^•^) 

I sinht I sinht j sinht I sinht j j 

= i(z^ - 1 ) [py{z)d.pp'’(z) - p-%)d,py{z)). 

For the minimally coupled massless case = 1 we have: 

= z(z 2 - l)Pf (z)P"*P(z) X 
2p 

~ |r[i + zp]|' 

2 sinh(7rp) 

TT 

(A.9) 


1 ip 1 ip 1 
- + -+ 


1 ip 1 
- + 


+ 


ip 1 


z — ip 21 + z 21 — 2 : z + ip 21 + 2 21 — 2 
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In fact, for v' ^ 'i this normalization is also valid. In [20] it is shown that one can expand 
the mode functions in the t —?■ 0 regime: 


1 

sinht 


(cosh t) 


2ip 

r[i — ip] 


^ip-i 


(A.IO) 


Using this expansion in (A.8) also results into the normalization (A.9). This normalization 
is valid for any t by the properties of the Klein-Gordon normalization. 


Bunch—Davies modes 

The Bunch-Davies modes are given in terms of linear combinations of the hyperbolic modes 
in (A.6). Schematically we have (using the orthogonality of the hyperbolic mode functions): 


- (X(T,P) Xo-,p) 


q=L^R q'=L,R 

,,h,R V ) 

NIp {o^pAq - l^p,qKq) ’ 

q=L,R 


(A.ll) 


where we used: 


(PP’^ )kg = )kg = 

(^pP.q^pp,q')^G = {P^^\P’^’^')KG =0. 

Using (A.7) we find: 

( a -U \ 8sinh7rp(cosh7rp- o-costtP) 

]T[u' + ip + lf' 

So finally we can substitute (A. 13) into (A.ll): 


— Npp 'Y^ {pp,q^p,q Pp,qPp,q) 
q=L,R 

2sinh7rp 8sinh7rp(cosh7rp — ucosttP) 
TT ^ |r[P + fp + l]|^ 

16 sinh^ Tvp (cosh Trp — a cos vrP) 

TT |r[P + ip + ijp 


(A.12) 


(A.13) 


(A.14) 


This is consistent with [20] , but note that we included an extra factor of 2 sinh np into the 
normalization, in order to simplify the expressions (A. 7). The normalized mode functions 
are the same as in [20], of course. 
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A.3 Mode functions for the massless scalar field 

Since we are mostly concerned with the massless minimally coupled scalar field [v' = 1), we 
state the normalized mode functions for that case explicitly in hyperbolic time coordinate 
t and in conformal time 77 = In tanh |: 


1 H 

Npp sinht 


P^^(cosht) 


H 

yj2p{p^ + 1 ) 

H 

\/2p(p2 + 1) 



{p csch t + i coth t) 


e {p sinh p — i cosh p ), 


(A.15) 


where we have chosen a convenient phase factor in the normalization, that does not affect 
the physics. The conformal time p is defined as: 


ds^ = ^ (y—dt^ + sinh^ t [dr^ + sinh^ rdVL^) 
f dt 

^ p = 


sinht 


, , t 

= in tanh -. 

2 


(A.16) 


Other useful relations between “hyperbolic” time t and “conformal” time p are: 

sinht =- - —, cosh t = —coth r/. (A.17) 

smhr/ 

At early times r/ —?• — 00 the mode function for the massless minimally coupled scalar field 
(A.15) behaves like a positive energy mode function: 


\/2p(p2 + i)' 


e ^' (psinhr/— z coshr/) 




\y/2p[p- + ij 

H{p + i) 

ismhr?— , =e 

\/2p(p2 + 1) 


(A.18) 


The infinite boost limit 7 — 7-00 (2.7) corresponds to the large t (or small p) and large 
momentum limit. In particular in terms of conformal time, the rescaling for small p reads 
p — 7 - pe~'^, implying that the combination pp is invariant in the limit. This gives 


e sinh ?7 — z cosh 77 ) = psinhue ( 1 -h 0 ( 77 ) | 

\ PV J (A.19) 

~ e-iPV ^prj - i ). 

This is exactly the mode function for a massless scalar field in the flat de Sitter slicing (up 
to the appropriate normalization). 
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B Power spectra for the massless field 


B.l Direct calculation 


Power spectrum in hyperbolic vacuum 

The power spectrum in the hyperbolic vacuum can be computed in a straightforward way 


^ ’ Im I'm' ^ 


X (Ohi (Vmdi’ + ikf’;"’) (Viw-p? ) isih) 


plm 


¥ , t\ 




yy ^pimXpi 


Sinh2(t) NpvNp,,, 

' ' Im /' yji' 


Im I'm' 
t2 


pipp-ip 

Vm|P //)=0 

PP' 


= Hp-p)- 




= Kp-p')- 


sinh^(t) 


yy i^zmi 




Im 
2 ^2 


_/V 2 

iVpp 


using 


bplm) b^p'l'm' — )bll/Smm' 




sinh2(t)27r2 iV2 


pp 


(B.l) 

where we used the completeness relation for the Ypim in the last step. For the massless 
minimally coupled case p' = 1 we have; 


(Pp|(/>p(/>p|Pp) = 6{p-p) 


P cosh^(t)+p^ 
sinh^(t) 47r2 (p2 -g l) 


(B.2) 


and for large t —)• oo 


(Pp|(/>p(/>p|Pp) 


5{p-p') 


p 

47r2 (p 2 + 1 ) 


(B.3) 


The power spectrum for the massless minimally coupled scalar field is given by: 


{Qh\4‘‘^\^h) 





dp' {^}h\4>p4>p\^h) 



47r2 p 2 _j_ 


p 2 

47r2 p 2 _|_ 1 ■ 


(B.4) 
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Power spectrum in Bunch—Davies vacuum 

The computation is similar to the previous case: 


Iml'm' aa' ,<y' 


X {^Bd\ i^c7plmXp,(j ^fjpim)(-P,<^ (^^a'p'l'm'Xp',a' ^a'p'l'm'^P'\^B d) 


i /2 


'^plm^p'l'm' JY ^ ’ {^BD\acrplm^]j'pri/m'\^BD) 

Iml'm' aa' Xp',^' 

H2 


T 


sinh2 t 
= ^ip-p') 


sinh2 t 


Ei^^'TE 


Im 


Xp,c 


N. 


Xp,cr 


using [Ufjpifn, ^a'p'l'm' 


— ^aa'^ll'^mm'^iP P ) 
,, ii 2 p 2 


= 5{p-p) 


sinh2 t 27r2 


E 


^ ( '^p,l'^p,L _j_ ^p,L^p,L 


^P,<T 


N 

CT 

. I 


N- 


py 


T' 

^pa 

I ^p,lPp,l pip pip I p-ip p-ip 

“I Ar 2 “I -i,,/ 


AT^ 


(B.5) 

In the last step we used the completeness relation for Ypim and the expansion of x in terms 
of the associated Legendre polynomials (A.6,A.7). Here we will compute the spectrum for 
the massless scalar field {u' = 1). For the massless minimally coupled scalar {P = 1) the 
cross terms involving P'Ppv and p-wp-w vanish: 


‘^P,lI^P,L ^ (g 7 rp_|_^)(g 7 rp_|_Q-^ ^ 

N^pa ^ cosh Trp + a 




(B.6) 


and similarly for the term involving P So for the massless minimally coupled case 

(z/' = 1) we have: 


/n ij. in \ w P^ I OippOipP + t^ppPpp 

(nsolMpliiBD) = s(p - p E -- 




P^p 


fj 2 2 

= Hp-p)zrp:2iPZ9Yl 


vr {e^P + (t)2 + (e '^p + a)'^ ^ip 2 


sinh2 1 27r2 16 sinh2 vrp cosh vrp + cr 


pW 


, 7/2 p 2 7rcosh(7rp) 

^ ^ sinh^ i 27r2 2 sinh^(7rp) 

i /2 p(cosh 2 (t)+p 2 ) 

= 0{p-p) _ o, , —TTTT-— cothfTrp). 


sinh 2 (i) 47r2(p2_|_i) 


(B.7) 


For large i —>■ oo we have: 


{^bd\4’p4>p'\^bd) = 5(p — p ) 


, //2 pcoth(7rp) 


4^2 p2 _j_ 


(B.8) 
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and 


{Q.bd\(I)^\^bd) = dp dp' {^Ibd\(I>p4>p'\^bd) 


dvr^ 


■ pcothvrp 

dp —^-h supercurvature modes 

p^ + 1 


(B.9) 


= / dlnp 


p^ coth TTP 


The power spectrum is given by 


^\bd{p) — 


p^ + I 

p^ coth vrp 


+ supercurvature modes. 


(B.IO) 


dvr^ p^ + 1 

which reduces for p 3> 1 to an approximately scale invariant spectrum: 

Aj(p) g. (B.ll) 

B.2 Reduced density matrix calculation 

In this section we derive the power spectrum in the Bunch-Davies state using an alternative 
method. We consider the reduced density matrix that remains after having traced out the 
degrees of freedom in the right hyperbolic patch. We find the same answer as in the direct 
calculation (B.10,B.7). The reduced density matrix for the left hyperbolic patch has been 
calculated by Maldacena and Pimentel [21] and is given by: 

PL,p,l,m = {\^Bd){^Bd\} 


= (1 - l7p)P) X] 


(B.12) 


n=0 


where for the massless scalar field 7 p and |n;j?,/,m) are given by^^: 

7 p(m = 0) = 


\n-,p, I, m) = 


plm' 

\/n! 




(B.13) 


For a point in the left hyperbolic wedge x € L the two point function is given by (2.22): 


{^bd\4‘p4>p'\^bd) — Tr ^^ {(Pp4>p'Pl} 


= d{p-p') — 




plm\ 
Im 


' sinh^(t) 

x(1-|t7)EI>7(2»+1) 


(B.14) 


= d{p-p')- 




p^ cosh^(t) + 1 + | 7 p 


sinh^(t)27r2 2p{p^ + l) 1 - 

r/ /4 P cosh2(t)+p2 

= o{j)-p) . / 2 , .X coth(7rp) 


sinh^(t)d7r2 (p^ + 1) 
which is equal to the result of the direct calculation (B.7). 


^^For the massive scalar field Maldacena and Pimentel apply a Bogoliubov transformation on the set of 
bpim operators to bring pL in the form of (B.12). 
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C Divergence of the energy momentum tensor 


As is the case in the Fulling-Rindler vacuum, the energy momentum tensor diverges at the 
null boundary of the hyperbolic patch. One could construct lightcone coordinates u = r] — r 
and V = 7] + r in order to calculate T^u- Equivalently, we consider the leading divergence 
of Ttt in the t —?• 0 limit, which is more convenient. 

Ttt = (C.l) 

For the massless case we have: 

{Tu) = + g^\decl>f + (C.2) 

One can calculate this directly using the hyperbolic mode functions (2.15) and the den¬ 
sity matrix (B.12) for the Bunch-Davies expectation value {Tfj^u)BD- Equivalently, for the 
leading order term we can use the Wightman functions G^{x, x') as given in [20]: 

{{Mf) = lim5i5f/G+(t,t'), (C.3) 

and similarly for the other coordinates. Note that the contribution of the supercurvature 
modes to the Wightman function only leads to subleading divergences^^. The contribution 
of the subcurvature modes to the Wightman function for the massless v' = 1 case is given 
by [20]: 


G+(t,t',C) = 


i/2 


1 


sinh t sinh t' Stt^ 


^ sinpC (cosh i-b ip) (cosh t'— ip) ^ tanh |- 


sinh Q sinh vrp 


l+p2 


I tanh I 


C = cosh r cosh r' — sinh r sinh r' (cos 6 cos 9' + sin 6 sin 9' cos{(j) — 4>')) . 

One can check the following: 

172 rPF 1 / 1 \ 

t', Q = ^ P^P^ + coth(7rp)^ -b O 

{{drctPf) = lim drdr'G^{t,t,Q = v-w / dp p{p^ + 1) coth(7rp)-^ -b O (t°) 

r'^r Jq ^ ' 

/"Pf 1 

{{decj))"^) = lim d 0 de'G'^{t, t, () = sinh^ r—^ / dp p(p^ -b 1) coth(7rp)-^ -b O (i°) 

O, ,0 X / 


(C.5) 


9’^9 


/■PF 1 

= \im^d(i,d^iG^{t,t,C) = sin^ 0sinh^ j dp p(p^-b 1) coth(7rp)^-b O (i°) . 

(C.6) 


'0 


^For u' > 0 the supercurvature mode contribution to the Wightman function is [20]: 

47r3 


Gt{t,t'X) = + l\V[v' -b (sinht sinht') 


/\ 1 ^' — I 


(C.4) 


For the minimally coupled massless case f = | the supercurvature mode becomes time-independent. The 
contribution to the energy momentum tensor is of subleading order. 
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Note that all these are divergent as t —)• 0, but they also show the usual UV-divergence. 
The UV-divergence is regulated by a cutoff pp- The difference {Ttt)H — {Ttt)BD will be UV- 
finite. We combine the components (C.6) to obtain {Ttt)BD- The expectation value {Ttt)H 
is obtained by replacing cothvrp —)■ 1, where we use the expectation value {bpim^pim +1) in 
the two different states 


(blimbpim + I)sz) = cothvrp 
+ I)h = 1- 


(C.7) 


Finally, we calculate the difference {Ttt)H — {Ttt)BD- 

^4 poo / 1 \ 

{Ttt)H - {Ttt)BD = ^ y P^P^ + 1) (1 - COthTTp) ^ + O j 


11 1 n 

2407r V 


(C.8) 


Note that we took the cutoff pi? to infinity and obtain a UV-finite integral. 


D The hyperbolic vacuum embedded in the Bunch—Davies state 


From [21] we have for the massless case v' = I: 



1 <8) \^h,r) 

(D.l) 

or suppressing the indices p,l,m: 





(D.2) 


with 7 = The left hyperbolic vacuum \^h,l) is not a state of the full system; we 

need information about the state in the right hyperbolic patch as well. The simplest way 
to embed the left hyperbolic vacuum in the full Hilbert space, we can consider the simple 
and symmetric state \^h,l) ® jUn,/?)- This state is not the natural vacuum state (the 
Bunch-Davies state) for the full de Sitter space. 

Proposition: 

\^h,l) ® \^h,r) oc “-|Us£)). (D.3) 

Proof: 

We will show that the right hand side of (D.3) vanishes when we act with any of the bpip^ 
annihilation operators. We use the expression for the hyperbolic annihilation operator bupm 
in terms of the creation and annihilation for Bunch-Davies modes (2.31), suppressing from 
now on the labels p, I, m: 

bl = ^ {oi'lcLa + • (D.4) 

a 

We want to show that bp acting on the RHS of (D.3) vanishes: 

E ^ I^bd) = 0. (D.5) 

a 

can calculate -I-1) by using the density matrix (B.12) 
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Consider the annihilation operator acting on (D.3): 


a±e-l>l“+“-|QBD) = [a± , |Obd) 


a± , |7p|a+ol 


=-|7p|a+al I 




(D.6) 


= -|7p|oLe-l>l“+“-|OBB). 


Substituting this result in (D.5) gives: 


E ^ (“l“" + |!2 bd> = Afp E 


Afv 


"I" ^ 


0^7 




iV- 


X-<T 


iVv 


"I" 


(D.7) 

It is easy to check that the quantity between brackets on the RHS of (D.7) vanishes for 
both cr = ±1. This finalizes the proof: 


bLe-\'r\->^-\nBD)=0 yp,l,m. 


(D.8) 


The state (D.3) is pure. Note that the symmetric and antisymmetric modes corre¬ 
sponding to cr = ±1 are entangled with each other and their reduced density matrices are 
thermal. 
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